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1 16. Abstract 'r =rarri mv-MCARF) for attached flow 

The multi-component sheet separation-flow 

is modified to accept the ptmaxI The viscous effects are in I 

gram (Analytical «®'^‘'oas, boundary layer dia- 

corporated into the o^riate source distribution. The 

placement thickness J ^ated into MCARF has been applied to 

separation flow model pressure distributions for 

single component airfoils. in close agreement with experi 
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SUMMARY 

The multi-component *^^^°^^fPg°^5^tex^hUt^separation 
flow is modified to accept the free vorte , effects are 

flow model of the ^^[^j^^tion by representing the boundary 

l;rer«cirnt"5hi^^^^^^^ aX apUriate source distri- 

bution. 

..e separation «o« -fl incorporated into ^ 

applied to single ti the stall are in close 

tributions for angles Even at higher angles 

agreement with experimental . trends of separation, de- 

lUltTt rirtlo^rtifieiti^a^rfncnise in pitching .cent 
coefficients are predicted. 

Although the program is ^“^^ned to hanole multl^^^ 

airfoils with checked out. It is suggested 

component airfoils such as the method of Nash and 

Sickfbe ln?c?lraSed into the model as opposed to the presen 

method of Truckenbrodt . 
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1.0 INTRODUCTION 


1.1 Background 

During the final design phase of an airfoil, it is 
essential to include the viscous as well as separation ef- 
fects at high angles of attack for the accurate evaluation of 
the aerodynamic performance. No general, mathematically 
closed form solution presently exists which describes the 
viscous separated flows around an airfoil. The availability 
of high-speed large capacity computers and the recent advance- 
ments in the potential flow and boundary layer computations 
has enabled the development of several computer programs 
which are based on iterative procedures between the viscous/ 
potential flow solutions for attached ai' veil as separated 
flows . 

In Reference 1, a computerized analytical model which 
computes the performance characteristics of multi-element 
airfoils in subsonic, viscous flow has been developed under 
a NASA contract to the Lockheed-Georgia Company. The model 
computes the viscous pres'^ure distributions, lifts, moments and 
local boundary layer properties on each element of an arbi- 
trarily arranged slotted airfoil in attached flow.^ The final 
viscous solution is obtained by an iterative technique by 
successively combining the potential flow solution with 
boundary layer displacement thickness. The details of the 
model and the computer program (MCARF) are presented in the 
next section. 

Multi-component airfoils operating at high angles of 
attack will usually experience regions of trailing-edge sepa- 
ration on the upper surfaces of one or more components. At 
low angles of attack, the cove regions of one or more com- 
ponents may contain a closed separation bubble which can re- 
duce the effectiveness of the downstream slot flow. In either 
situation, the performance of the system is reduced from that 
which would be expected with completely attached flow. In 
practice, some amount of separation is to be expected, and 
there is a need, therefore, for the capability to predict the 
performance of multi-component systems having regions of sepa- 
rated flow. 

Until recently, our ability to model separation was 
limited to models that simulated separation by source distri- 
butions employing the assumption that the pressure everywhere 
in the separated region is constant (Ref. 2) . None of the 
methods calculate the pressue in the separation region directly , 
but rely on some criterion to determine the separated flow 


pressure level. In general, the methods predict the upstream 
pressure distributions in a satisfactory manner once a suit- 
able source "outflow*' has been chosen. 

More recent approaches (Refs. 3 and 4) employ an inverse 
method to determine the effective airfoil shape to achieve 
constant pressure in the separated flow region. There is 
some question as to the uniqueness of the solutions obtained 
by this procedure. Since none of the methods use a physically 
realistic wake model, it is difficult to see their applica- 
bility in the stall region where a direct analysis procedure 
for the prediction of aerodynamic characteristics of airfoils 
is desired. 

A procedure which does employ a realistic wake model for 
separated flow is described in Reference 5. This method, 
developed by Analytical Methods personnel, uses free vortex 
sheets to separate the free stream fluid from the wake region 
(a region of lower total pressure) . The strength of each vor- 
tex sheet and its location are determined as part of the over- 
all calculation. The procedure allows for a direct calculation 
of pressures in the separated flow region. The details of the 
analysis and the computer program (CLMAX) are presented in the 
next section. 


1.2 Present Approach 

The primary objective of the present work is to modify 
the MCARF program to accept the free vortex sheet separation 
model developed for the CLMAX program. The viscous effects 
are incorporated into the calculation by representing the 
boundary layer displacement thickness with an appropriate 
source distribution over the airfoil surface. This approach 
enabled incorporation of the separation wake model into MCARF 
with a minor modification to MCARF, since it does not change 
the panel geometry of the airfoil surface. For the present 
study, the program is capable of analyzing a multi-component 
configuration with separation only on one of the components, 
presumably the trailing component. The effect of the sepa- 
rated region will, of course, impact the flow around the 
other components . 

The inclusion of the separation model into MCARF involves 
the development of two iteration loops. The inner (wake) itera- 
tion loop for the potential flow relaxes the wake shape for a 
prescribed separation position, while the outer (viscous/poten- 
tial flow) iteration loop predicts the separation point loca- 
tion and the displacement thickness distribution from the 
boundary layer analysis for each of the components. 
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2.0 NOMENCLATURE 


A 

C 


m 


sep 

IPOT 

ITR 

ITRW 

ITRWMX 

K 

N 

s 

V 

WL 

WF 

x,z 

a 

Y 

P 

AH 


Aerodynsunic influence coefficient 
Chord length 
Pressure coefficient 

Lift coefficient 

Moment coefficient about the quarter-chord axis 

The value of the boundary layer shape factor at whxch 
the flow is assumed to separate 

Choice of potential flow index 
Viscous/potential flow iteration index 
Wake shape (inner) iteration index 
Maximum number of wake shape iterations 
Kernel function 
Number of panels 

Coordinate along airfoil surface 

Velocity 

Wake length 

Wake fineness ratio 

Cartesian coordinates 

Angle of attack 

Stream function 

Vorticity strength 

Density 

Increase in total pressure in the separation region 
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3.0 DESCRIPTION OF COMPUTATIONAL MODELS AND PROGRAMS 

3.1 Multi-Component Airfoil (MCARF) Program 

In this section, a brief description of a computerized 
analytical model and the computer iiiti-element 

predicts the performance characteristics of the multi ele^^^ 
airfoil in attached, subsonic viscous flow are presented. T 
model computes the viscous Pressure distributions, lift and 

moments, Ld boundary layer properties ^J^ouf so?u- 

ayV>^<-rarilv arranoed slotted airfoil. The final viscou., »uxu 
tion is obtained by an iterative ^®<=hnique which ^mbines an^ 

inviscid solution with the boundary ^ layer closed 

naoa 'Phs» surface of each element is approximated as a ciosea 

polygon with segments represented by c5^o|\athe-” 

larities. The boundary layer solution is comprised of roatne 

matical models representing s^®te-of-the-art technology 
laminar, transition, and turbulent boundary 

tional boundary layer model has been incorporated to compute 
the characteristics of a confluent boundary layer which 
fleets the merging of the upper-surface boundary layer witn 

the slot efflux. 


3,1,1 Program Operation and Capa bilities 

The airfoil program is composed of three 
(1) geometry specification; (2) potential ^ 

arl layer. After data input and geometry specification, ^e 
program enters an iterative cycle ”hich involves the deteraina 
tion of the interrelationship between the potential flow and 
tho honndarv laver After each iteration, a convergence check 
is made which insists of a simple comparison of the computed 
ontQ Experience has shown five iterations are 
lirt a converged solution. This rapid convergence 

necessary to obtain ite?ative technique in which the physical 

can be ^ each iteration is modified taking 

ISrISoSL boundary layer thickness computed from the previous 
iteration. 

more 

IS 

SSfUSSircISuStirStltrSShSSSrpSrsure^are 
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potential Flow Sol -ularity distribution method 

liJlFrSeniHtS^^^^^^ 

equation 


° , X _ o z(s) sin (ct) 

= U x(s) cos (a) ^00 

00 


( 3 . 1 ) 


• ♦- r(s C) is the 

«... . u r%r.»rrJs’S€S?- ;'“L? 




Equation (3.2) becomes 


4 , - 



U. 


cos 


(a) - Zi sin 


(a) 




(3.2) 


A. . is the velocity induced at 

r;r.e": 

“S^"%S%rx!o:LrrnSlSoIltion formula as used. 


- . ''1-1 I°-I ~ °i-l l_!_Ili °1*^ ' 


(3.3) 




Sj +1 ' ®j-l 


For this method, “ ®“^''^°'J?he®MhSns'’L2®N-f nSb^ of 

^Vs'-aS *. ^^'^rVeclafJta condi^ T tof riose spac- 

r„f iel“Se ?ramnred.e. This %t« -"dition 

^S!rfs that the J?f Se4r 'Ih4 upper\nd lower surface 

|5\ir«fiUnrSof anl”Sat at the trailrn, ed,e 


(T,) “ 


- (Y^), 


(3.4) 


•4-v, +-v,in new Kutta condition 

,oovirercar,fhf sriSr-^^^ for%Ssp-ed,ed as well as 

supercritical airfoils. 

Boundary layer solution 

The pressure coefficients account for com- 

SSlsiS 
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ard input to the boundary layer portion o£ th^^^^ 

boundary layer consists of an layer is 

confluent boundary layer. The turbulent regions. The 

?? airis oSe of the inigue features of this program. 

Ks; ;y:: 

!J?Sf“SfSIed\S^?h“;itr^ Of CoAe^ 

friLfer^ii^t!’'5oiunrBLS 

frsf-ioj* 

An initial check is made to determine ir criterion 

layer is stable or unstable based . boundary 

eskblished by SohUctrng 

pKfafly Snivel 'transition prediction curve 
?lfrccu?red, the inibiel quantrtres needed. to^^start the tur 

bulent calculations * ^r^na bubble with correspond- 

curred, the with corresponding reattach- 
ing laminar stall or a short transition loca- 

After computing the transition location and the 

of the mixing between the slot to the trail- 

forward element, and can the pressure distribu- 

ing edge of “e element, depending upo^the^^^ Goradia (Ref. 7) 

If S;sef In'Sf rslSm3 Lf SI Je^ 

tfSiSf fnllflf ln“fiafrnallgous to that for a free Jet 
flow. 
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Bcuivalent UrfoU Geometry ^ 

me vISloS^IStitM Sr.Sif steps are as toXXows. 

n 4 *ial flow solution for the basic 

(1) Compute a potential 

airfoil. properties based on the poten- 

(2) Compute boundary lay P 
tial flow solution. 

“ rrr,T::::r“:r:. ». 

“?lomanS%oefflcients is obtained. 

in the present version ^Jn'’m?p™icSnrwith'’the\ 

-ifriSl 

thickness ^^uivalent airfoil geometry, 

the computation of equivalent 

3.2 cT.MAX Program 

In this section a method °^ig“desSibed! ^It is^an 

airfoils up to and beyon^ ^o^ential ^ pcteSial 

sf,uire In inner 

Uera?ioi^?o%s?alish their shapes. 


. a. it... r-,f +-hp Flow Model 
3.2.1 nascript ion or >-ug 

^ flow’ field with separatiov, is shown in Figure 1. 
reaions are identifiable. 


Several 


n • «n 1 • The Potential Flow Region 
to the fluid. 
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_ _ I __ 5 , The Boundary Layer 

Region 2. me owu j v,«a hiqh shear, 

The thin ^®^^°’^^*'®ggea°which ^^Jfate Tignlflcan? vi?tifit?!' 
and hence, viscous stresses wnx 


Region 3: The Free Shear Layer 


3 n 3: Tne tree laver 

The thin flow region TO?ti?ity“«<>"?P«‘ 

signifies*'^ • 


Region 4: 

The 
with low 


The Wake 


wake between the two 

vorticity and insigni 


shed boundary 
ficaiit viscous 


layers is 
stresses • 


a 


region 


T^ARic Assu mptions 

c *-vio flow is defined by these 
An approximate model of th 

assumptions. 

(1) The boundary layer ^5®®hence, can be represented 

significant th ohness an^ ^rrl^iinea acroaa whioh 

it veiocity. 

,2) The wake doea not have =ig;i«i|snt^vortrcity^an^^ 

P°-hant tota 


The flow field car he ~rrifcftie^ St^Sfared^^f a°“ 
;U?“iS°d!rtPbufion^Sith a ?^if’“hld“ihrmwjust Se 

boundary layer ana 

The mathematical bouSdarf coSdUioS ^fmet. 

?jr^fs^itrof of ?he'ro 4 ?irtrshret\epreaentin^ the free shear 

layer is not known a priori. 
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Boundary Conditions 


The boundary condition for 


the airfoil surface is flow 


tangency or 


^ • n = 0 


(3.5) 


«here S = unit nomnl vector, and ^ is the total vorticity 
vector. 

When allowing the >=oundary layer displacement effect, the 
right side of this equation will be non zero. 

The free vorticity sheets are located on streamlines and 
there is no static pressure drop across them. 


up py-ovi mat ions for the Free Shear Layer 
(1) Wake Shape 

initially, the ® rit«~L“y,'ltS 

of free shear layers "‘"st e calculations in which the 

from an initial r wy iteration suggested the 

vortex sheet shapes were obtained by iteration gg 

initial shape shown i" ^"^'^^?h^uppe?ail SSef Sheets are 

basic ^^!^P^°?frcurveI pKsL? from the separation 

represented by parabolic ^ ^ slope at the up- 

points common to a point . free stream direction and 

Stream end is the *«ean between t^ free, stre^^ 

I'll.lTtl ?SI*re?n"iSriine! dlatanc^ WL downstream from the 
wake midpoint (Figure 2) . 


(2) Wake Length 

Early calculations indicated that jesults^were^sensi-^^ 

tive to the of the fr ^ith relatively 

with experimental ^esul 0.1c to 0.2c beyond the 

short wakes; i.e., wakes extena ng length, WL, 

Ja?e^bSnin5;stigped for a iJ^wa^^^ 

i'S!! SS"!s oKa'Iy mStiplying the "height" of the wahe 
(Figure 2) by the wake fineness ratio, WF. 



FIGURE 2. INITIAL WAKE GEOMETRY 


(3) Wake Pressure 

The approximation /express io^^ total 

free shear layer is /e? the strength of the free vor- 

KrlSStt? «5lrS/iru|pe. Shear layer, if the average 
velocity in the layer is denoted by 

/ \ (3.6) 


V 


'“I 


^outer ■*■ ^inner) 


then 


V ^ = V + Y„/2, and 

^outer u 


V . = V - Y../2 f 

^ inner u 


. . . _ V - V . , on the upper sheet . 

since the vorticity, Yu “ ''outer inner _ _ 

(The vorticity in the lower shear layer is Y^^ - ^^ner " outer* 
The jump in total pressure across the shear layer is t en 


= ”inner " “outer 


= 1 Pinner " ) ‘ 1 

-I Pouter 


^ 1 
. J 


= - pVYu = 


(3.7) 


... 

same everywhere. 


( j 
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3.2.2 Potential Flow Solution 


The boundary condition of flow tangency on the airfoil 
surface gives the integral equation: 


^■sy(s)ds+ (f Kde - f Kde ) 

T. U 


+ V. 


n = Vj, (3.8) 


Where the constant value 

niSil JelJoit? component 0.1 ^Int 'is- " 

MOlltSd Slt^the’ eluent , 6s, ^^^®/tne"ai?fSl°Sid“he 

till ^d llltk I:- S-il^r-vrcnS'cn"’'^ 

reeSM°prl^Sgd= 1 and VL. ,-|PecUv-r laSir 
former is a function of position on the airroii, 
is a constant. 

iSSsirnSrar-iooftrind-^ 

distribution. 

Numerical Solution 

. - 4 .„„v- -ic rporesented by an inscribed polygon, 

Figurri.“‘5hrLS!;!dSal panels representing the polygon each 

ha?e a linear variation of y°-^“^i^gj='f;aneis “'unOor^i^ 
tex sheets are represented by a ntmber p 

"°""f^s"denMy"vi in ^qn. (3.9) 

panel noted y unknown sequence, {yi>f 

?h2 lnte,?al equations in ‘^e unknown function, theSfare 

Slrunknow^"?fva5Lr(S? STnels) , iSt'the auxiliary condi- 
tions remove two unknowns: 

at the upper surface separation point, Yggp ~ 

(y^ being the value at the lower separation point) , 



Also, the Y value just ?SirarlVl^S^TOvm 

^“n?S?cing ?L eurS« boundary condition (^n. <3.9 )) 

are amenable to solution. 


Calculation of Pressures 

Having found the vorticity, the velocity at gy J;" 

the flow field can be evaluated by adding to the free stre^ 

Se velocities iSuced by the vorticity and source distributions , 
we p?ess«S are oaloulLed from the velocities^aocordrng to 
the Bernoulli equation which is expressed non-dimensionally as 


(v«) 


(3.10) 


P - Pc 


^here C = / q„ = = increase in total pres- 

sure ovL that”at infinity. Note that AH = 0 everywhere except 
in the wake_region for which it was previously shown (Eqn. (3.8)) 

that AH = pVYy • 


3 2.3 Structure of the Iterative Pr ocedure 

Figure 4 shows an outline of the method which has two 
iterative loops. 

(1) Wa)ce Shape Iteration 

ticity IS ass^ea CO divided into a number of 
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INPUT 


ATTACHED 

POTENTIAL 

FLOW 

SOLUTION 


BOUNDARY LAYER 

analysis 

(SEPARATION POINTS 
& BOUNDARY LAYEIi 
SOURCE DISTRIBUTION) 


V/AKE A/IODELING 


POTENTIAL 
FLOW SOLUTION 
(WITH WAKE 
VORTEX SHEETS) 


ITRW < ITRWAK? 


C£ CHANGE 
< 1 % 

OR ITR > 8? 


MORE INCIDENCE 
CASES 


(POTENTIAL FLOW 
PRESSURE 
1 DISTRIBUTION) 


FIGURE 4 . 


STOP 


PROGRAM OUTLINE 
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The wake shape is calculated as follows. Using the previous 
vorticity distribution, velocities are calculated at the panel 
midpoints on the free vortex sheets. The computed velocity dis- 
tribution is then used to arrive at the new wake shape by re- 
orienting the wake panels in such a way they are tangential to 
the local flow directions, starting from the • separation points. 

The upper and lower sheet downstream end points, which 
were coincident in the initial wake, are allowed to move in- 
dependently in subsequent iterations. At each iteration, the 
wake influence coefficients at the surface control points are 
recalculated, and a new potential flow solution is obtained. 


(2) Viscous/Potential Flow Iteration 

This outer iteration loop takes the potential flow pres- 
sure distribution over to the boundary layer analysis and re- 
turns with the separation points and with the boundary layer 
source distribution. The source distribution is determined 
directly from the boundary layer solution as 

a = £ (0^«*) 

where Ug is the streamwise potential flow velocity at the edge 
of the Boundary layer, and 5* is the displacement thickness. 

The addition of this source distribution modifies the normal 
velocity, Vjj (Eqn. (3.8)), at each panel control point. The 
sources are set to zero in the separated region. 

The program generates a new wake shape using the -new 
separation points together with information from the previous 
iterated wake. A new potential flow solution is then obtained, 
and so on. The outer iteration is terminated when the change 
in Cj, is below 1%. A limit of twelve iterations is currently 
imposed within the program. 
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4.0 MODIFICATION OF MCARF WITH SEPARATION PLOW MODEL (MCARFM) 

The primary °|i.Sax fJS'Sortef 3heet°se|f«tlon 

SeL'’'?r?Mreirt!?l a'aescrlphlOh of this modlficatioh rs 
presented. 

For the present study, the am is capable^of^analyzi^^^ 

a multi-component =??^ 9 uration with s P t^sk of the present 

components, the “ailing component^^^The^ ^3^ iteration loops to 
investigation involves the a® P MCARF. The inner (wake) 
incorporate the separation model in MCARF. i ^ 

iteration loop for the “Surthfluter (viscous/ 

a prescribed separation ts the separation point 

potential flow) ^°2mpnt^thickness distribution from the 

location and the < 3 xspl®cement thick ^^j^ponents. The concept 

boundary layer analysis Z- reoresent the effect of the 

of displacement thickness is ^^ed to represent 
SLious viscous layers on the outer Potential 

the VISCOUS effects of t following two options, 

potential flow oalculation^ y displacement thickness is 

In the first option (IPOT ^ arrive at an equivalent 

combined with the Ji® ous%otential flow iteration, 

airfoil geometry after each_viscous/p^t^ displace- 

In the second option " aeometry, a distribution of sources 

ment thickness to Jil i 2 ed^fot the simulation of the 

along the airfoil surface is the so-called surface 

viscous flow ‘3isplaoe;;ent effects^ framework of thin boundary 

transpiration method which, witnin method of adding 

layer theory, is oo^P^ete y q to the basic airfoil 

geometrically the ‘displacement thickness 

is ^S^vlntlnf t^rdlpJ^Ihe s4psration flow modal into 

MCARF . 

4 . 1 Development of Separation Flow Mode J. 

The major code changes Jj^e'^^POTVL^the^pot^ 

involve the modification (Over lav 2) For the present 

flow solution »^o;^tine in to either a single corn- 

investigation, ’^‘^®.'"°f®J;„lti-component airfoil with separation 
ponent airfoil or to ^ nnlv The separation 
restricted to the trailing component) • is predicted 

location over the upper surpce ( (Hn „ith the pressure 

from the boundary "jj^jy^^^^tenfiarflow U or 

distribution computed from iteration. It is assumed 

separated) solution in th P surface always occurs at the 

that the separation at the lower sutta 

trailing-edge of the airfoil. 
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initially, tha sttaamllnas are He'^preSn^d 

SjSlaSS!!^ ^u?4r;LIS l?£ e II a on Poln» 

“Mill tre"fnrltr£a U=n ana^tte^local surlace_^ 
slope. The common point ^®'^aa“*am P miapoint (Figure 

waki line, alstance WL lie heighi of the 

like S%lflaMeilsf I 

11‘inruf Thi lillliliinllnla^lSte aeveloping the 

CLMAX program. (See Appendix) 


4 . 2 Potential Flow Solu tion 

For the "g- JjqulKitJ^distr iLtio^me^^ Oellers. 

is obtaineo using the singular y details of the method 

The governing flow equation J 3 . 2 ^,^and^^the^detai^^^ 

were presented in Section * vortex strength is assumed 

into N sa^ments over each of which J°rtex^streng^ 

to be a constant. The free ^n*~c of uniform 

■4;,Krr»H.sr 

this source sheet. 

in MCARF, the panel i2i°"”'^iirnol|ltla°il"llltinl^AlRir"’' 
bounaaries ana J^il the lame panel Information is 

iliriSoughout lie complete visoous/potential flow iterative 

scheme for an attached flow the upper surface 

flow case, the separation point ”^ 3 ^®atbitrl?y parameter 

vortex sheet leaves ^he. air foil surface another 

which varies from one viscou^potential flow ire 

sSfef '°?h?r?eprne?tJrscheme f^a temporary arrangement at 
lach of the visc?us/potIntial flow iterations, and the changed 
panel geometry is replaced with the original panel geome y 
before proceeding to the next iteration. 
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For the separation , th®the*sipatation%S^ con- 

(E,n. (3.M) o£W is replacea^wrth^the sep^^,^^ ,„,aratl=aUy 

dition, Y-ep = » ana uoper surface separa- 

for the last four segment ‘^°^J®g®trailing edge of the airfoil. 

Son point and the She slparation region is 

Hence, the unlfom vortrcrty ‘^sep> equations consisting 

not an additional unkno^^ in t Y Additionally » m the 

o£ N+1 unknown values of y e ^ ^ downstream of the =?P“»J“Jre 
separation regj^on,^ values veiling edge of the airforl 

point and at the upper ^ 

Assumed to be zero (see Figure jj. . . 

in the present j°5f°reSS^etwSJn’"ShS‘^StJeSrf^^ 

satisfied by ‘^adiacent segment midpoints, J^gpiption 

e^imtSaSL^thrunSnown value of^^ SSS^^N^boSary conditions 

eauations with N unknown values or t ^ similar com- 

lirfoil surface. For the separated flow | exceptions. 

Stational procedure is adopted with the „ake region 

Th^kSta condition is Jon/^?hl vorticity values on 

as discussed earlier this separation point and the 

the downstream segment ^ are assumed to be zero, 

trailing-edge segment on the '^PF eliminated from the original 

Snci Swo dlditional equations are^^slrmin^^^ H -2 equations 

set of N equatxons to arrx As discussed earlier / ^ 

“f; rs^rr-aariifona? u^lknown ani 

p!Ii?! \hf rreat^f lis^rlbririn ?S separation region is 

computed using Eqn. (3.10). 


4.3 c;^rncture of the _Xj ierative Sche me 


4 , 3 brrucx, Lii-^ 

The iterative scheme adop^d ^J^^^^/^FigSre'^lH^with a 
very similar to that of th^^^^ it was assumed that a 

few exceptions. the iteration when changes less 

solution is converged at the iteration; however, the 

than 1% from the value at the limited to 8 due to 

maximum number of limitations in the P^°^^®^'it can 

the storage and data ^rans input by the user, and ^^.can 

the present program, criterion demands, m the original 

be as high as any convergence criteri additional 

ScABF program, f aJsSmsd to be converged when 

restriction that the s ^ i n the previous iteration, 

is within 0.005 from “®siSg!e’'airmulti-element airfoil cases 
“S«^1rc°o"nvrrredi°lbb!ins after as few as 5 iterations. 
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,.3.1 «aS^shaES_(in!HEL^S^^ 

The Iterative loop for ^ton^S„i J . «iihi^^ls l°°P4^e 

involves the P°tentral flo«_^=°i^„t, which was computed^^^^^_ 

the upper „tion of the are represented 

f„rtfaSy^th”^%ra^'roSn^^ 

bv parabol-c curves this section. ^ 'i^nels (NSEPWK) • 

ita-;r«gf;;a;;.r.a » S.K. 
sara%r,;s~^»s;.:”aafi^^ 

rr?^"a., ®|P«k\s a |inplo inaeP-de^ is a 

^ ^nQT_T_V bs A-i--ar-*V 

lunction of the anple o a mop is Pre^nted^in^^^ 

Pi,ur^e ?rn»e\ 

previous vorticity free vortex she . starting 3*^ 

fotnrsr«Srw|£c°inoiae^t^r^^^^^^^ 

obtained. frn’RWMX) is an input 

The inaximum “f JssSs at f “|ff„,lxim« of three 

leparation on s a =°"terged a £ixed 

f p«sin?s i«tf =t/r=%^5?!°lnras thi pSent 

rhire“i-o|«gea afte| |?^^ 3 ,To«-rr. the iser can change r 
program » liKwmA 
if needed. 
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FIGURE 5. INNER (WAKE) ITERATION LOOP 






criterion of Nash and ^ ^fairly ^ simple” to^a^^ routine 

routine TURB and it ^„ith^the separation flow model. 

Lto the present MCARF program ^^^|^^tion criterion of 

For the present value for HSEP is arrived at 

r.reraSrniertir^^^^^^^^ -- • 

In the surface transpiration ^method ^^°gjterthe^displace- 

i.:nreflect°Sf ?hfviscous' layers The strength a of this 

^qSLfle" source is calculated from 


a = 


A/ 


P 

e 


A 


(4.1) 


1 •+...« Is the coordinate along 

where p is density, displacement thickness, and e 

the airfoil surface, 6* is tne ^ . e respectively. The 

» refer to local and edge of the boundary layer, 

local conditions are at the |ion modifies the normal 

in the separation region. 

The outer iteration scheme starts wrth^ 

Tn fn-fal^^Kr^fpr^ . 

distribution, using this ^nformatio followed to ODtain 

is defined and the inner ^teration p solution with 

Ihe relaxed wake shape and the potential separated 

separated wake. The boundary y location and the boundary 

yields a new ^^^Kl vtscSSs/^otential flow iterative 

layer source distribution. gj-gent solution is obtained, 

process is repeated until a con | point and lift 

Fiaure 7 shows the history of the P . ^ attack of 19.06 . 

coSficient of a GA(«)-1 “jJ^^^efce'ohiJIoteriatics of the present 
This demonstrates a good conve g 
model . 


4.4 


^ the present investigation, a major 
During the course the presen gg^eral viscous/ 

difficulty was f^'^^^^tered. While separation 
potential flow Jteration , continued to move towards the 

location over the upper sur co velocity distri- 

leading edge "'°"°tonically. A cl ^ the shape factor, 

S“‘rcrrarefl^t^rr%fpaS ^ 

rt!Sn“o?"a/ilfl?ia;s‘forwIra of the separation location fro™ 
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Iteration 




Iteration 




the previous iteration, where the separated 

was first thought that this might be due to the pesent Kutta 
condition, whi?h assumes that the vortex strengths 
ticallv for the last four segment corners near the upper sur 
iJIatiL aSd the. lower surface. of the trailing^edge and ^f^the^^ 
vortex strength in the separation wake^ ^sepu ^sep£^ 

strength at the lower surface trailing edge. . f 
alternate Kutta conditions such as linear vortex 

strength for the last two segment “^„=1^^r^‘p~'t1on 

Strength at the last of the trailing edge were tried, 

location and the __ of steep rise in velocity 

These conditions ^ion somewhat, but did not elim- 

gradients near the separat investigation, it was con- 

inate it completely. Based increase of velocity gradients 

eluded that this problem of model and the 

may be due to the P°^®JJ^^^Jj.j^Jlation of the^aerodynamic matrix 
procedure adopted in tn flow boundary conditions 

In MCARF. .in of equa?Lns between pairs 

are satisfied by formulating a thereby eliminating the unknown 

of two adjacent segment "'J;‘^P°^ oe While the model proved 

stream function value over the surface, wnii 

to be very effective for t^e slpara?ed flow and 

rmori"dei!ilerand\Eo?ou^^ investigation is needed to conclude 
whether it is indeed the problem. 

the velocity distribution at the last two segme 

separation region such mint ahead of the 

are equal to the 9^^dient at the third seg location 

separation region. ®JPPfJ^^ Ihe vLcous/potential flow 

to move forward well . ^ consistent way. Admittedly, 

iterations progress in a ^ purely scientific basis, 

if f u gof enfinL^ing Approximation, 
axpeoiafly lippof ed by the experience gained while developing 
the AMI CLMAX program. 

In the MCARF program, HSEP, the value of the 
shape^?acf r at .wh?ch . the flow separates was jssum^ ^1.8. 

Sf®a??ackrSIf ifesAigated varying the value of HSEP 
angles of ettack) were /,-oncludcd that the optimum value 

fif 5 ISp = l! St toteS on the correlation between the computed 

and test results. 
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AS discussed earlier J°th2”viscSus^^e?S!^^n^ 

two options boundary layer displacement thickness 

first option (IPOT = 0) , the boundary layer aip 

ifused'^to arrive at an equivalent J^^^S^^IJSSI-option (IPOT = 1) , 
viscous /potential flow iteration. appropriate source 

the viscous effects are simulated by adding^an^ 

distribution along the otpion is used. However, 

tion flow modific tion, the a^ion, it was discovered 

during the course of the second option was 

that MCABF uses ^he first option ly a ^ investigated, 

not completely checked out. In fa ' substantially larger 

it was found that the second option the first 

Cp values for attached tlow , ^gj-ijnent. This discrepancy 

edition correlated well with th P following reasons: (i) 

may be attributed to one or both of the toiio 

the maximum wake centerlines in the presence 

inadequate ^®P5®^?^ihution This problem does not present any 
of the source ‘^^^tribution. s boundary layer calcula- 

difficulty for separated . location and the source 

tion is terminated at the separat separation region. Hence 

strengths are assumed to angles of attack (with 

it is suggested to use option 1 gr^low ang ^ separated 

no separation) , es are made in the present program 

flows. The appropriate option. 

(MCARFM) to incorporate this particular 

5.0 DISCUSSION OF RESULTS 

The method was ^PP^j®^ 3 ^^°irfoS^shape"'?ep?^ a difficult 

NACA 4412 airfoils. b^tions are available from experi- 

test case and Pressure distributions are ^^tack. 

ments at NASA-Langley for a range or angie 

The first set of results. Figures 8-11, ^are^f or ^ 

potential flow iterrtLns. 

coefficient for this cas , Figure 9 compares the cal- 

a good convergence * distributions at a = 20.05°, 

culated and ®^P®tin'e'^tal P again, there is good 

which is 3 "st beyond the St a^l.^^0 and experimental 

agreement. Figure 10 P _ . ,.o^ The comparison is not as 

pressure distributions predicted separation point 

good as the other cases at about 0.1c. The rate 

has not reached the ®^P® Calculated separation point with 
of forward movement of the ^ , movement of the separation 

21.14- is slow. The sudden fo^ari^movemeut^ 

SfattiSk^lompKeTiSures 9 and 10) is difficult to predict 
in this case. 28 
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figure 8. 


r. 1 ? r>jvT rULATED AND EXPERIMENTAL PRESSURE 



^ Experiment 

——Present Calculations 


COHP«.ISO« OF CM-COLMED A® EXFEFI«E«Ttt PBESSORE 

DISTOIBOTI0N„0N^A GAWHI , 0 . 15 , 


1 





COMPARISON OF CALCULATED A^ _EXPERIMENTAI- PRESSURE 

°e i “1 Ic'loS and H - 0.15) 





a, degrees 


O Experiment 

Present Calculations 


a, degrees 
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a, degrees 


O Experiment 

Present Calculations 


a, degrees 



FIGURE 12 , 






a, degrees 


COMPARISON OF CALCULATED AND EXPERIMENTAL 
LIFT CHARACTERISTICS FOR A NACA 4412 AIRFOIL 
(Re = 6.3 X 10®) 


FIGURE 13. 
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Figures 11 and 12 show the lift and J 

of the GA(W)-1 airfoil at Reynolds ® es they 

X 106 , respectively. As can be seen from these ^J^ures, they 
show good agreement with experiment. Figure shows the lift 

coefficient comparison with eperiment ' f ^_Jee 

4412 airfoil at Re = 6.3 x lO^. The calculated values agree 
very closely with the experimental values up to 

in all of the cases- studied under the present 
the cLculated values compare well with the results of the 

CLMAX program. 


6.0 CONCLUSIONS 

The separation flow model incorporated into MCARF has been 
applied to single “tfoils^ M^thfstall‘’ltrin close 

f^?reieS^wJ?h S^rLe^tal -asurements Even at Wgher angles^^ 

fr, tJ^?1oe&fts,nn“in«ersrin pitching .cent coefficients 
are predicted. 

aitfc^ifSiS fe;a?::!c a^%^rira!?iS 

Hicks be incorporated into the model as opposed to the prese 
method of Truckenbrodt. 
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APPENDIX; PROGRAM PAR AMETERS 

In the present separation the'^dlscussion 

parameters is the the SMI CLMSX Pt°9tam, 

frfnS S? airiili sections and thicjcness ratios wer^ ir^re^ented 
fn%»5!^ ?^if d°a?a“if SsId^affSI^^rnce guide for the 
present investigation. 

The input paraiaeters and the £°™f *^'^Stth’'SrLle?tiM 
(MCSETM) is identical with the input of program 

°lv«!af Sk °?Se iS^t JUilwes and the format specification 
on CARD A are presented below. 


CARD A; 


Coliomn 


6-10 


11-20 

21-30 


Separation Flow Variables 

Variable Description 

NSEPWK Number of panels in the separa 
tion region 

ITRWMX Maximvim number of wake shape 

iterations 

WF Wake fineness ratio 

Hcpp The value of boundary layer 

shape factor at which the flow 
is assumed to separate 


Format 


FIO.O 

FIO.O 


wake Length = VJF x Wake Height (see Figure 2) 



1 / (Thickness/Chord Ratio) 


FIGURE 14. 


DEPENDENCE OF WAKE FINENESS RATIO, WF, 
aTRPOTL THICKNESS /CHORD RATIO. 
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